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We report a fully inclusive measurement of the flavor changing neutral current decay B ^ Xs'y 
in the energy range 1.7 GeV < i?^ ™ '^ < 2.8 GeV, covering 97% of the total spectrum, where c.m.s 
is the center of mass system. Using 605 fb~^ of data, we obtain in the rest frame of the _B-meson 
B{B Xsj : > 1.7 GeV) = (3.31 ± 0.19 ± 0.37 ± 0.01) x 10"", where the errors are statistical, 
systematic and from the boost correction needed to transform from the rest frame of the T(45') 
(c.m.s) to that of the B-meson, respectively. We also measure the first and second moments of the 
photon energy spectrum as functions of various energy thresholds, which extend down to 1.7 GeV. 
The results are preliminary. 

PACS numbers: 12.39.Hg, 13.20.He, 13.40.Hq, 14.40.Nd, 14.65. Fy 



Radiative _B-meson decays may offer a view of phenom- 
ena beyond the Standard Model of particle physics (SM) . 
In the SM, these decays proceed via a flavor changing 
neutral current (FCNC) decay, which consists of a loop 
process. Yet to be discovered particles, such as charged 
Higgs or supersymmetric particles, may be produced vir- 
tually in the loop and produce a measureable deviation 
from the branching fraction predicted by the SM. 

The predictions of the branching fraction at order 
(NNLO - next to next to leadiiig order) (3.15 ± 0.23) x 
10-4 [H, (2.98 ± 0.26) X 10-4 ||| and the average of ex- 
periment measured values (3.55 ± 0.26) x 10-^ are in 
tacit agreement. An updated experimental measurement 
would further test this agreement, and, moreover, give 
stronger constraints on extensions to the SM e.g. Min- 
imal Supersymmetric Standard Model Q and left-right 
symmetric models Q . 

The photon energy spectrum is also of great impor- 
tance. At the parton level, the photon is monochromatic 
with energy E « mb/2 in the 5-quark rest frame. The 
energy is smeared by the motion of the &-quark inside the 
B meson and gluon emission. A measurement of the mo- 
ments of this spectrum allows for a determination of the 
&-quark mass and of its Fermi motion. This information 
can then be used in the extraction of the CKM matrix 
elements \Vcb\ and \Vub\ from inclusive semileptonic B 
decays Q . A measurement of the low-energy tail of the 
photon spectrum is important in this context Q. 

Belle has previously measured the B ^ Xg'y branch- 
ing fraction with 5.8 fb-^ and 140 fb-^ of data using 
semi- inclusive ^ and fully inclusive approaches p^ . 



respectively. Other measurements include those from 
CLEO [HI and BaBar [ll,[il,[i3. 

Here we present an update of our fully inclusive mea- 
surement [10|, based on a much larger dataset and with 
significant refinements, which includes an unfolding of 
detector effects on the measured spectrum that improve 
the measurements of the branching fraction and spectral 
moments, respectively. We also extend the photon en- 
ergy range to i^^ ™ '^ > 1.7 GeV, covering more of the 
spectrum than ever before, where c.m.s refers to the cen- 
tre of mass system, which is equivalent to the rest frame 
of the T{AS). 

The B Xgj decay is studied using the Belle detector 
at the KEKB asymmetric e^e~ storage ring The 
data consists a of sample of 604.6 fb^^ taken at the T(4S') 
resonance corresponding to (656.7 ±8.9) x 10^ BB pairs. 
Another 68.3 fb-^ sample has been taken at an energy 
60 MeV below the resonance and is used to measure the 
non-BB background. Throughout this manuscript, we 
refer to these data samples as the ON and OFF samples, 
respectively. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer described in detail elsewhere [l^. The main 
component relevant for this analysis is the electromag- 
netic calorimeter (ECL) made of 16.2 radiation lengths 
long CsI(Tl) crystals. The photon energy resolution is 
about 2% for the energy range relevant in this analysis. 

The strategy to extract the signal B Xgj spectrum 
is to collect all high-energy photons, vetoing those orig- 
inating from 7r° and 7] decays to two photons. The con- 
tribution from continuum e+e^ ^ QQ il — u, d, s, c) and 
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QED type events is subtracted using the OFF sample. 
The remaining backgrounds from BB events are sub- 
tracted using Monte Carlo (MC) distributions scaled by 
data control samples. 

Photon candidates are selected from ECL clusters of 
5x5 crystals in the barrel region (—0.35 < cos 9 < 0.70, 
where is the polar angle with respect to the beam axis, 
subtended from the direction opposite the positron beam. 
They are required to have an energy E'^^ ™ *' larger than 
1.4 GeV. We require 95% of the energy to be deposited 
in the central 3x3 crystal array and use isolation cuts to 
veto photons from bremsstrahlung and interaction with 
matter. The center of the cluster has to be displaced 
from any other ECL cluster with E > 20 MeV by at least 
30 cm at the surface of the calorimeter, and from any 
reconstructed track by 30 cm, or by 50 cm for tracks with 
a measured momentum above 1 GeV/c. Moreover, the 
angle between the photon and the highest energy lepton 
in the event has to be larger than 0.3 radians at the 
interaction point. 

In the Belle detector, a non-negligible background (1%) 
is due to the overlap of a hadronic event with energy 
deposits left in the calorimeter by previous QED in- 
teractions (mainly Bhabha scattering). Such compos- 
ite events are completely removed using timing informa- 
tion for calorimeter clusters associated with the candi- 
date photons. The cluster timing information is stored 
in the raw data, and is available in the reduced format 
used for analysis only for data processed after the sum- 
mer of 2004. This divides our data set into 253.7 fb~^ and 
350.9 fb^^ samples of reprocessed data without and with 
timing information, respectively. To minimise composite 
background due to Bhabha scattering and two-photon 
processes that contaminate both T(4S') and continuum 
data samples, we veto any candidate that contains an 
ECL cluster with energy exceeding 1 GeV within a cone 
of 0.2 radians in the direction opposite our photon can- 
didate as measured in the c.m.s frame. In the second 
data set only photons that are in time with the rest of 
the event are retained. The efficiency of this selection 
on signal events is larger than 99.5%. We veto candidate 
photons from tt'^ and rj decays to two photons by com- 
bining each B Xsj candidate photon with all other 
photons in the event. We reject the photon candidate if 
the likelihood of being a tt*^ or 77 is larger than 0.1 and 0.2, 
respectively, these yield, on average, background suppres- 
sion factors of 4 and 2, respectively. These likelihoods are 
determined from MC and are functions of the laboratory 
energy of the other photon, its polar angle 9 and the mass 
of the two-photon system. 

In order to reduce the contribution from continuum 
events, we use two Fisher discriminants calculated in the 
c.m.s frame. The first discriminant exploits the topol- 
ogy of B —^ Xgj events and combines three energy flows 
around the photon axis. These energy flow variables are 
obtained using all particles, except for the photon candi- 
date, we measure the energy in the three regions defined 
by e < 30°, 30° < e < 140°, 6 > 140°, where 9 is the 



angle of the particle to the candidate photon. The second 
exploits the spherical shape of BB events and is built us- 
ing ten event-shape variables including Fox- Wolfram mo- 
ments [17'] for the full event and for the partial event with 
the photon removed, the full- and partial-event thrusts 
and the angles of the thrust axis with respect to the 
beam and the photon direction. To optimise these se- 
lection criteria, we use a MC simulation 18] containing 
large samples of BB, qq and signal weighted according 
to the luminosities of the ON and OFF samples. In the 
optimisation step the signal MC used is generated as in- 
clusive B Xsj and exclusive B K*j. The inclusive 
component Xg is defined as a resonance of spin-1 with a 
Breit-Wigner form and a mass of 2.4 GeV/c^ and width 
1.5 GeV/c^. The X, system is hadronised by JETSET 
and subsequently re weighted to match the prediction of 
the DGE model [H with mb(MS) = 4.20GeV/c2, 
with the mass extending no lower than 1.18GeV/c^ to 
agree with the corresponding world average branching 
fractions 3. To improve the understanding of the pho- 
ton energy spectrum at low energies, the selection criteria 
are optimised to maximize the sensitivity to the signal in 
the energy bin 1.8 GeV < E^-'"-'' < 1.9 GeV. 

After these selection criteria we observe 4.15 x 10^ and 
0.25 X 10^ photon candidates in the ON and OFF data 
samples, respectively. The spectrum measured in OFF 
data is scaled by luminosity to the expected number of 
non-BB events in ON data and subtracted. The formula 
used to subtract continuum background is as follows: 

^SBj^^c.m.s(ON)-j _ ^ONj^^c.m.s(ON)-j 

ON gON 

^Hadronic B^Xsf p 

" ' ,OFF ■ ,OFF ' 

''Hadronic '^B^X.f 

•iVOPF(Fsi?---^(OPF)) 

where e is the efficiency of Belle's hadronic selection [23| 
or of this analysis' {B Xsj) selection criteria in con- 
tinuum events at either ON resonance {^/s = 10.58 GeV) 
or OFF resonance {^/s = 10.52 GeV) energies, and a is 
the ratio of ON to OFF resonance integrated luminosity 
corrected for the energy difference (a = 8.7557(±0.3%)). 
The factors Fe and Fn compensate for the slightly lower 
mean energy and multiplicity of particles in OFF com- 
pared to ON events. We find Fn = 1.0009 ± 0.0001, 

Fe = 1.0036 ± 0.0001, = 0.9986 ± 0.0001, and 

ON 'Hadronic 

'gr/°^ = 0.9871 ± 0.0014. The ON and scaled OFF 

spectra and their difference are shown in Fig. [1] 

We then subtract the backgrounds from B decays from 
the obtained spectrum. Six background categories are 
considered: (i) photons from tt*^ —f 77; (ii) photons from 
77 — > 77; (in) other real photons (mainly decays of lu, 77', 
and J/tpi E^nd bremsstrahlung, including the short dis- 
tance radiative correction (modelled with PHOTOS 21]); 
(iv) ECL clusters not due to single photons (mainly -ft'^'s 
and n's); (v) Electrons misidentified as photons and; (vi) 
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FIG. 1: ON data (open circle), scaled OFF data (open square) 
and continuum background subtracted (filled circle) photon 
energy spectra of candidates in the c.m.s frame. 



Contribution 


Fraction 


Signal 


0.190 


Decays of tt" 


0.474 


Decays of rj 


0.163 


Other secondary 7 


0.081 


Mis-IDed electrons 


0.061 


Mis-IDed hadrons 


0.017 


Beam background 


0.013 



TABLE I: Relative contributions of the BB backgrounds after 
selection in the 1.7 < £'^ ™ 7( GeV) < 2.8 range 



beam background. The spectra of the background of pho- 
tons from B-meson decays with respect to the expected 
signal is shown in Fig. [5] and listed in Table HI The net 
background of this type is a factor five greater than the 
signal. 

For each of these categories we take the predicted back- 
ground from MC and scale it according to measured 
yields wherever possible. The inclusive B —> ir^X and 
B i]X spectra are measured in data using pairs of pho- 
tons with well-balanced energies and applying the same 
ON— OFF subtraction procedure. The yields obtained 
in data are on average 10% larger and 5% lower for tt'~' 
and rj than MC expectations. The observed discrepancy 
between the measured and simulated tt" 77 spectra is at- 
tributed to the branching fraction assumptions used for 
the generator Beam background is measured using 
a sample of randomly triggered events and added to the 
BB MC. 

For each selection criterion and each background cat- 
egory we determine the K ™'^-dependent selection effi- 



71"— >yy 
T]— >yy 

Other decays 
Beam bkgd 
Mis-ID e 
Mis-ID hadron 
Signal 




FIG. 2: The spectra of photons from _B-meson decays passing 
selection criteria as predicted using a MC sample. 



ciency in OFF-subtracted ON data and MC using ap- 
propriate control samples. We then scale the MC back- 
ground sample according to the ratio of these efhciencies. 
The efficiencies of the tt" and 77 vetoes for photons not 
from tt" and rj are measured in data using one photon 
from a reconstructed tt" , where the other photon of the tt'~' 
is excluded from the search over the remaining photons 
for the next best 7r° or ?/ candidate (highest 7r° or 77 like- 
lihood) . Consequently the best formed tt" or i] candidate 
used in the calculation of the likelihoods is most likely a 
random combination, and therefore suited to measuring 
the effect of the vetoes. The tt" veto efficiency is mea- 
sured using a sample of photons coming from measured 
7r° decays. We use partially reconstructed D*~^ £)°7r+, 
Z?" K'TT^-K^ decays where the 7r° is replaced by the 
candidate photon in the reconstruction. The rj veto effi- 
ciency for photons from tt'^'s and event-shape criteria effi- 
ciencies are measured using a 7r° anti-veto sample, which 
is made of photons with a 7r° likelihood larger than 0.75 
(ie, no TT*' veto) and passing all other selection criteria. 
Other efficiencies are measured using the signal sample. 
Beam background is negligible after the application of 
the OFF time veto. In the sample of data where the 
veto is unavailable we scale the background according to 
a comparison of yields between MC and data for high 
energy (i?^ ™ '' > 2.8 GeV) photon candidates found in 
the endcaps of the ECL. This sample after continuum 
subtraction is a clean sample of ECL clusters from beam 
backgrounds. 

The ratios of data and MC efficiencies versus i*;^^ '"-^ 
are fitted using first or second order polynomials, which 
are used to scale the background MC. Most are found to 
be statistically compatible with unity. An example is the 
effect of the 7r° veto on photons from tt^s that escape the 
veto in the partially reconstructed D* sample, which is 
shown in Fig. [31 

An exception is the efficiency of the requirement that 
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FIG. 3: (LEFT) The tv° veto efficiency in the partially recon- 
structed D* sample for both Data (circles) and MC (squares) 
and (RIGHT) their ratio fitted with a first order polynomial. 



95% of the energy be deposited in the central nine cells 
of the 5x5 cluster, which is found to be poorly modelled 
by our MC for non-photon backgrounds. We estimate 
the efficiency for data using a sample of candidate pho- 
tons in OFF-subtracted ON data after subtracting the 
known contribution from real photons. This increases the 
yield of background (iv) by 50%. The yield from the six 
background categories, after having been properly scaled 
by the above described procedures, are subtracted from 
the OFF-subtracted spectrum. The result is shown in 
Fig. 21 After these subtractions the yield in the spec- 
trum above the endpoint of B decays is compatible with 
zero, 1245 ± 4349 candidates. 

To measure the branching fraction and the moments 
we correct the raw spectrum using a three step proce- 
dure: (i) divide by the efficiency of the selection criteria 
i.e. the probability of a photon candidate passing cuts 
given a cluster has been found in the ECL, as a function 
of the measured energy in the c.m.s frame; (ii) perform 
an unfolding procedure based on the Singular Value De- 
composition (SVD) algorithm [2^, which maps the spec- 
trum from measured energy to true energy thereby un- 
doing the distortion caused by the ECL; (in) divide by 
the efficiency of detection i. e. the probability that a pho- 
ton originating at the interaction point is reconstructed 
in the ECL, as a function of the true energy. Data are 
divided into 50 MeV wide bins. Step (ii), which was 
not performed in our previous analysis, is essential for a 
consistent extraction of partial branching fractions and 
moments as a function of lower energy thresholds. The 
unfolding matrix, derived from signal MC, is calibrated 
to data using the results of a study of radiative di-muon 
events, which gave the ECL response in data and MC in 
an energy and acceptance range consistent with our anal- 
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ysis. We use five sig nal models: KN 
DGE [21], BBU ^ and GG Values of the parame- 
ters of the signal model used in the signal MC are derived 
from fits to the signal spectrum shown in Fig [4] The two 
error bars for each point show the statistical and the total 
error, including the systematic error which is correlated 
among the points. In order to obtain the total B Xs"f 



branching fraction we apply corrections for the contribu- 
tion from Cabibbo suppressed B — Xdj decays. The 
ratio of the B Xgj and B — > Xdj branc hing frac- 
tions is assumed to be R^/s = (4.0 ± 0.4)% [30|. We 
apply corrections to derive the measurements in the B- 
meson rest frame, using a toy MC approach. We generate 
photon 4-momentum in the rest frame of the i?-meson 
using signal models referred to earlier, and generate B- 
meson 4-momentum using their known fixed energy and 
1 — cos 9^ distribution in the c.m.s. Repeating this ex- 
ercise many times yields photon energy spectra in the 
rest frame of the i3-meson and the c.m.s, from which we 
extract corrections used to yield measurements in the B- 
meson frame. The correction is derived as a mean over 
all signal models while the root-mean-square is assigned 
as the uncertainty. After correcting for the acceptance 
we derive distributions of the partial branching fraction, 
first moment (mean) and second central moment (vari- 
ance) of B Xsj as measured in the c.m.s and B rest 
frame for lower energy thresholds as shown in Fig. (5] In 
the range from 1.7 to 2.8 GeV in the rest frame of the B- 
meson, we obtain a partial branching fraction, and the 
first two moments of the energy spectrum: 

B{B ^ X,j) = (3.31 ± 0.19 ± 0.37 ± 0.01) x 10"^ 



2.281 ± 0.032 ± 0.053 ± 0.002 GeV 



(E^) - {E^f = 0.0396 ± 0.0156 ±0.0214± 0.0012 GeV^ 

where the errors are statistical, systematic and from the 
boost correction, respectively. 

The full results, the systematic error budget and corre- 
lation coefficients for five lower energy thresholds {Ei^ = 
1.7, 1.8, 1.9, 2.0, 2.1 GeV) are listed in Table|TTl The to- 
tal systematic error is derived from a sum in quadrature 
over all sources. We vary the number of BB, the ON to 
OFF ratio of integrated luminosities and the correction 
factors applied to the OFF data photon candidates and 
assign the observed variation as the systematic associ- 
ated with continuum subtraction. The parameters of the 
correction functions applied to the 7r° and rj yields are 
varied taking into account their correlations. As we do 
not measure the yields of photons from sources other than 
tt'^'s and rj's in BB events, we independently vary the ex- 
pected yields of these additional sources by ±20%. For 
the model dependence in correcting for the acceptance we 
use four signal models in addition to the default model, 
and assign the maximum deviation from the default as 
the uncertainty. The error on the photon detection effi- 
ciency in the ECL is measured to be 2% using radiative 
/i-pair events, and also affects the estimation of photons 
from BB. For the uncertainties related to the unfolding 
procedure, we vary the effective rank parameter up and 
down by one in the SVD algorithm. 

In conclusion, we have measured the branching frac- 
tion and photon energy spectrum of -B ^ Xgj in the 
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TABLE II: The measurements and correlation coefficients of the branching fraction, mean and variance of the photon energy 
spectrum for various lower energy thresolds, i?;^, as measured in the rest frame of the B-meson and the contributions to the 
systematic uncertainty. 



B{B ^ X,7) (10- 



(E^) (GeV) 



AE^ = (E?,) - {E-,Y (GeV") 



E^ [GcV] 



1.7 



1.8 



1.9 2.0 



2.1 



1.7 



1.8 



1.9 



2.0 



2.1 



1.7 



1.8 



1.9 



2.0 



2.1 



Value 
zbstatistical 
zbsystematic 
±boost 



3.31 3.24 

0.19 0.17 

0.37 0.24 

0.01 0.01 



3.12 2.94 2.62 

0.15 0.14 0.12 

0.16 0.12 0.10 

0.02 0.02 0.05 



2.281 2.290 2.305 2.326 2.350 

0.032 0.025 0.019 0.015 0.011 

0.053 0.028 0.014 0.007 0.005 

0.002 0.002 0.004 0.005 0.006 



0.0396 0.0350 

0.0156 0.0096 

0.0214 0.0081 

0.0012 0.0005 



0.0292 
0.0058 
0.0027 
0.0008 



0.0227 0.0170 

0.0033 0.0017 

0.0009 0.0006 

0.0009 0.0012 



Systematic Uncertainties 



Continuum 
Selection 

Other B 
Beam 
resolution 
Unfolding 
Model 
7 Detection 



0.18 
0.20 
0.07 
0.24 
0.02 
0.01 
0.01 
0.03 
0.03 
0.01 



0.11 
0.15 
0.05 
0.13 
0.02 
0.01 
0.00 
0.02 
0.02 
0.01 



0.08 
0.11 
0.04 
0.06 
0.01 
0.02 
0.00 
0.01 
0.01 
0.01 



0.07 

0.08 
0.02 
0.02 
0.01 
0.02 
0.01 
0.00 
0.00 
0.01 



0.07 
0.06 
0.01 
0.01 
0.01 
0.03 
0.01 
0.00 
0.00 
0.01 



0.030 
0.023 
0.012 
0.033 
0.001 
0.006 
0.002 
0.005 
0.005 
0.000 



0.016 
0.012 
0.006 
0.016 
0.001 
0.005 
0.001 
0.003 
0.003 
0.000 



0.008 
0.006 
0.003 
0.007 
0.000 
0.005 
0.001 
0.002 
0.002 
0.000 



0.004 
0.003 
0.002 
0.002 
0.000 
0.004 
0.001 
0.001 
0.001 
0.000 



002 
001 
001 
000 
000 
004 
002 
000 
000 
000 



0101 
0114 
0075 
0124 
0006 
0009 
0014 
0014 
0014 
0001 



0.0040 
0.0039 
0.0023 
0.0051 
0.0003 
0.0006 
0.0008 
0.0006 
0.0006 
0.0000 



0.0012 
0.0014 
0.0007 
0.0017 
0.0001 
0.0005 
0.0006 
0.0002 
0.0002 
0.0000 



0.0004 
0.0005 
0.0003 
0.0004 
0.0000 
0.0004 
0.0003 
0.0000 
0.0000 
0.0000 



0.0004 
0.0001 
0.0001 
0.0000 
0.0000 
0.0004 
0.0001 
0.0000 
0.0000 
0.0000 



Correlation coefficients (combined statistical and systematic) 



1.7 



1.8 



AB 
1.9 



2.0 



2.1 



1.7 



1.8 



{E-,) 
1.9 



2.0 



2.1 



AEi 



1.7 



1.8 



1.9 



2.0 



2.1 



1.7 
1.8 
AB 1.9 
2.0 
2.1 
1.7 
1.8 
1.9 
2.0 
2.1 
1.7 
1.8 
1.9 
2.0 
2.1 



1.000 0.959 
1.000 



0.811 
0.942 
1.000 



0.699 
0.839 
0.939 
1.000 



0.604 
0.720 
0.823 
0.959 
1.000 



0.455 
0.269 
0.031 
-0.004 
0.023 
1.000 



0.322 
0.129 
-0.107 
-0.143 
-0.107 
0.967 
1.000 



AEi 



0.114 
-0.073 
-0.291 
-0.332 
-0.296 
0.838 
0.946 
1.000 



-0.083 
-0.251 
-0.442 
-0.494 
-0.476 
0.636 
0.793 
0.942 
1.000 



142 
291 
464 
531 
548 
489 
645 
824 
954 
000 



0.848 
0.807 
0.680 
0.612 
0.546 
0.342 
0.155 
-0.066 
-0.230 
-0.252 
1.000 



0.857 
0.878 
0.817 
0.767 
0.689 
0.149 
-0.047 
-0.280 
-0.438 
-0.438 
0.945 
1.000 



0.722 
0.822 
0.869 
0.863 
0.795 
-0.094 
-0.290 
-0.516 
-0.660 
-0.642 
0.782 
0.935 
1.000 



0.528 
0.678 
0.814 
0.870 
0.848 
-0.252 
-0.431 
-0.640 
-0.779 
-0.777 
0.581 
0.782 
0.946 
1.000 



0.445 
0.568 
0.700 
0.846 
0.910 
-0.174 
-0.329 
-0.530 
-0.696 
-0.751 
0.497 
0.677 
0.840 
0.942 
1.000 



energy range 1.7 GeV < E'^-'^-'' < 2.8 GeV in a fully 
inclusive way. For the first time 97% of the spectrum 
is measured [3l| allowing the theoretical uncertainties 
to be reduced to a very low level. Using 605 fb~^ of 
data taken at the T(4S') and 68fb~^ taken below the 
resonance, we obtain B{B Xgj : > 1.7 GeV) = 
(3.31 ± 0.19 ± 0.37 ± 0.01) x lO""*, where the errors are 
statistical, systematic and due to the boost correction, 
respectively. This result is in agreement with the latest 
theoretical calculations The results can be used 

to place constraints on new physics [1^ and determine 
SM parameters such as the 6-quark mass [SS*! . 
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FIG. 4: The extracted photon energy spectrum of B ^ X^^d'y- 
The two error bars show the statistical and total errors. 
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FIG. 5: (TOP) Partial branching fractions, (MIDDLE) mean, 
and (BOTTOM) variance of B ^ Xs^y in the (LEFT) c.m.s 
and (RIGHT) and in the rest frame of the S-meson for lower 
energy thresholds. The two error bars show the statistical 
and total errors. 



